Antiproton and Radio Constraints on the Dark Matter Interpretation 
of the Fermi Gamma Ray Observations of the Galactic Center 
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Recently, it was suggested that the gamma rays observed by the Fermi Gamma Ray Space Tele- 
scope from the direction of the galactic center could surprisingly well be described by a dark matter 
annihilation scenario, both in terms of their angular distribution and energy spectrum. Here, I point 
out that such a scenario would be in considerable tension with existing antiproton data, in particu- 
lar the new PAMELA measurements. Radio data are even more constraining, thus disfavoring the 
dark matter hypothesis and making an astrophysical explanation of the observations much more 
plausible. 



The innermost region of the Milky Way is one of the 
classical targets for indirect dark matter (DM) searches 
In fact, the expected high DM density close to the 
galactic center would make it the single brightest DM 
annihilation source in the sky. This apparently favorable 
situation for the hunt for DM signals is, however, ob- 
scured by the fact that the region is also astrophysically 
very rich and complex; as a consequence, it will generally 
be difficult to unambiguously distentangle a DM signal 
from the resulting, still not very well understood back- 
ground 0. 

In the past, the potential observation of a DM signal 
from the galactic center region has already been claimed 
a couple of times - in particular in gamma rays Q , but 
also in microwaves [H or the annihilation radiation from 
positrons Q ■ While there still might remain a bit of con- 
troversy in some of these cases, evidence is certainly not 
compelling and the general picture has been that more 
refined analyses and new data tended to disfavor the DM 
hypotheses previously put forward @ (see, however, 0]). 
The most recent in this series of claims is the 'possible 
evidence for dark matter annihilation' seen by Goodc- 
nough & Hooper They find that the gamma rays 
from this direction, as observed by the Fermi Gamma 
Ray Space Telescope @, can well be described by a DM 
annihilation scenario where the DM particle has a mass of 
x ss 28 GeV and annihilates purely into bb final states, 
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with a cross section of (av) « 9 • 10 _26 cm 3 s _1 . For 
the dark matter distribution in the halo, a profile with 
p x ( r ) k r -11 has to be invoked, i.e. somewhat steeper 
than the standard NFW profile [10(. Rather than argu- 
ing for or against the plausibility of this scenario, and 
whether or not potential astrophysical backgrounds were 
adequately taken into account, let us here check the in- 
ternal consistency of this proposal. 

Apart from gamma rays, antiprotons are also an in- 
teresting means of indirect DM detection as they are 
inevitably produced whenever it is kinematically possi- 
ble and the annihilation products have color charge 11 1 . 



Unlike gamma-rays, however, they do not travel unper- 
turbed through the galaxy but are deflected by inhomo- 



FIG. 1: This Figure shows the antiproton flux as measured by 
BESS '95/'97 Qj], CAPRICE QJ], BESS polar QJ] and, most 
recently, PAMELA [l7| . The yellow band shows the predicted 
flux in secondary antiprotons, taking into account the full 
range of uncertainty in the propagation parameters []]§] . The 
dotted line indicates the primary antiprotons produced in the 
DM model by Goodenough & Hooper Q] and the thick red line 
shows the resulting total flux as calculated with DarkSUSY 
[Isl ]. assuming the smallest allowed secondary flux. 



geneities in the galactic magnetic fields. While this intro- 
duces a certain amount of uncertainty, the main aspects 
of antiproton propagation are well understood and can 
nicely be described in terms of semi-analytic diffusion 
models where the parameters are adjusted to fit existing 
cosmic ray data, in particular the observed boron over 
carbon ratio B/C [l2|. As a result, the expected back- 
ground flux of secondary antiprotons, mainly produced 
in the collisions of cosmic ray protons with interstellar 
matter, is tightly constrained, even well beyond the cur- 
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FIG. 2: The shaded region indicates, roughly, the combina- 
tion of DM mass m x and annihilation cross section (av) that 
could in principle describe the gamma ray spectrum observed 
by Fermi from the direction of the galactic center; the cross 
corresponds to the best fit model obtained in [i§] which was 
also used in Fig. [T] Models above the solid line are excluded 
due to the overproduction of antiprotons; the dashed (dash- 
dotted) line shows the constraints from radio observations, 
assuming a cutoff in the DM density profile at r c = 1 pc 
(rc < l(T 2 pc). 



rently accessible energies. For illustration, Fig. [T] shows 
the expectation for this background flux, fully taking into 
account the uncertainty in the propa gati on parameters, 
as predicted roughly three years ago [ijj]. The remark- 
able agreement with the subsequent data of BESS polar 
[III and (still preliminary) PAMELA ^ provides an 
important test for the underlying diffusion model. 

In the same Figure, the total antiproton flux in the DM 
model proposed by Goodenough & Hooper is indicated, 
as calculated by using DarkSUSY [l8| with its default set 
of propagation parameters (implementing the procedure 
described in [19|). As it stands, the resulting antiproton 
excess at energies below about 5 GeV would clearly not 
be compatible with the data. One has to keep in mind, 
however, that the B/C analysis of the allowed range of 
propagation parameters is limited to sources in the disk; 
most of the annihilations, on the other hand, take place 
in the halo, so one cannot expect the predictions for the 
flux of primary antiprotons (produced directly by DM 
annihilations) to be as tightly constrained as for secon- 
daries. By adopting extreme sets of propagation param- 
eters, still compatible with B/C, the flux of primary an- 
tiprotons could, in fact, be up to an order of magnitude 
smaller (or larger) than what is shown in Fig. [Tl 12|. 



Including the new (PAMELA, in particular) data sets 
in the analysis, on the other hand, would constrain the 
allowed range of propagation parameters even more _2(| . 
For a given annihilation cross section, the annihilation 
signal enhancement that is expected due to the effect of 
DM clumps in the halo (often referred to as "boost fac- 
tor") is, furthermore, probably stronger for antiprotons 
than for gamma rays, thus making the antiproton bounds 
somewhat more constraining than discussed so far. The 
reason for this is that for a spiky DM profile like the 
one considered here, the largest contribution to the total 
flux in gamma rays comes from the innermost part of the 
galaxy, where clumps most likely have been completely 
disrupted by the strong tidal forces; antiprotons probe a 
much larger portion of the halo and the total flux is not 
very sensitive to the region close the center. Taking into 
account these considerations, I will in the following re- 
quire that the primary flux from DM annihilations alone 
should not exceed five times the maximally expected flux 
of background (i.e. secondary) antiprotons; while a more 
detailed analysis of the uncertainties involved would cer- 
tainly be interesting, it is clearly beyond the scope of the 
present letter. The resulting constraint in the (av) vs. 
m x plane is shown in Fig. [2] and falls well below the re- 
gion that could describe the Fermi gamma ray data in a 
DM setup like the one proposed in [8[ . 

By changing the propagation model, or adopting ex- 
treme propagation parameters that are barely excluded 
by the B/C analysis, one could probably still argue for 
the possibility to evade the antiproton constraints pre- 
sented in Fig. [2] However, even more powerful con- 
straints can be derived by considering electrons and 
positrons as DM annihilation products and the synchro- 
ton radiation they emit when propagating through the 
galactic magnetic field. For the galactic center, it has 
earlier been pointed out that this usually provides more 
stringent bounds than for gamma rays, or photons at 
other wavelengths, especially for low mass DM models 
21 1 . In the monochromatic approximation for the syn- 
chroton emission the resulting luminosity from a cone 
towards the galactic center can be written as f22| 



uL(v) = 2ir 



(av) 



drr z p x 2 (r)E e± N(E e ±). (1) 



In this expression, E e ± = y/ Aitm z e v /B = 
0.43 GeV(^/GHz)5 (B/mG)"5 is the energy of an 
emitted electron or positron at which synchroton 
radiation of frequency v is emitted and N(E) is the 
number of electrons and positrons, per annihilation, 
above that energy. The strongest constraints now derive 
from comparing the above luminosity to radio data at 
the lowest observed frequency, v = 0.408 GHz [23| (the 
observed region corresponds to a distance of ~ 0.14pc 
from the galactic center). 

Before being able to derive these bounds, one needs 
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of course also a model for the magnetic field in the in- 
nermost galaxy. As an upper bound on its strength, 
one may consider it to be in equipartition up to the 
accretion radius R acc « 0.04 pc of the black hole, 
with magnetic pressure and kinetic energy of the in- 
falling matter balancing each other, which leads to 

B(r < i?acc) « 7.2 mG (0.04 pc/r) 5 / 4 [H; above this, 
conservation of the magnetic flux dictates a scaling 
B(r) (x r~ 2 , which in the following will be assumed down 
to the typical, constant galactic value of B ~ ^G. As the 
opposite extreme of a rather small magnetic field, one 
may assume it to stay constant below i? acc . 

Following thus the same procedure as in [i^], the re- 
sulting constraints on the DM mass and annihilation 
cross section are also included in Fig.O Here, an artificial 
cutoff of the DM profile was introduced at r c = 1 pc, be- 
low which the DM density was assumed to stay constant. 
This is an extremely conservative assumption since such 
a cutoff would usually only be expected to be set by the 
annihilation rate 25j , at a scale several orders of magni- 
tude smaller than this. As also indicated in the Figure, 
a smaller cutoff would further strengthen the radio con- 
straints by a factor of up to about 200 (for r c < 10~ 2 pc, 
these constraints actually do no longer depend on the 
cutoff for the halo profile considered here). 

It should also be stressed that the radio constraints 
shown in Fig. [2] are only very weakly dependent on the 
assumptions about the magnetic field - in contrast to the 
assumptions about the DM profile which enters quadrati- 
cally in Eq. ([T]) - since the main contribution comes from 
not too small radial distances, where the uncertainty in 
the magnetic field is not very big; varying the magnetic 
field configuration between the two extreme cases men- 
tioned before leads to barely visible differences in the 
bounds. Note also, from Eq. J!}, that the luminosity 
actually increases for smaller values of B, so the (some- 
what counterintuitive) assumption of a more a less con- 
stant magnetic field at the /^G level extending to the very 
central part of the galaxy would lead to even stronger 
bounds. 

At this point, one may wonder whether there is a 
chance to modify the proposed characteristics of the DM 
model in such a way as to change the general conclusions 
we have reached so far and that are nicely summarized in 
Fig. O a DM interpretation of the Fermi data is in strong 
tension with existing antiproton and, in particular, radio 
data. While the expected flux in antiprotons is not very 
sensitive to the form of the DM profile, for example, the 
radio constraints in principle arc. However, the profile 
cannot be changed significantly because of the need to fit 
also the angular dependence of the Fermi data down to 
an angular resolution of about 0.1° (which corresponds to 
a distance of ~ 10 pc from the galactic center; note that 
a cutoff at this scale would be completely ad hoc and not 
have any astrophysical motivation already the cutoff at 
lpc included in Fig. [5] is artificially large and only added 



for comparison). The only point left to consider, then, is 
the possibility of modifying the annihilation channel. For 
kinematic reasons, however, the annihilation to Higgs or 
gauge bosons is not possible; annihilation to other light 
quarks produces essentially the same spectrum in p, e 
and 7 as for bb, and gamma rays radiated from light lep- 
ton final states would have a very different energy distri- 
bution compared to what would be required in order to 
reproduce the Fermi data. This leaves only one poten- 
tially interesting channel, i.e. annihilation into r^, but 
one may easily check that also in this case the resulting 
spectrum in gamma rays is too hard to fit the data in a 
satisfying way. 

To summarize, the recent antiproton data very nicely 
confirm earlier expectations about the flux of secondaries 
and thereby provide interesting constraints on exotic con- 
tributions like from DM annihilation (a fact that has 
been already observed, at higher energies, in the ear- 
lier released p/p data [26j]); especially at low DM masses 
these constraints on the annihilation cross section arc 
quite relevant as they are on the order of what is ex- 
pected for thermally produced WIMPs. Radio data are 
usually even more constraining. Taken together, these 
constraints are in strong tension with the possibility that 
the gamma rays observed by Fermi could be due to DM 
annihilation in the galactic center. The radio bounds are, 
in fact, so stringent that they seem to seriously obstacle 
even the discrimination of a potential sub-dominant DM 
contribution to the signal. 
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